1 



Update on 

C. Bernard? * T. Blum? T. DeGrand? C. DeTar? Steven Gottlieb? U. M. Heller* J. Hetrickf 
C. McNeile, d K. Rummukainen, A. Soni, b R. Sugar? D. Toussaint 8 and M. Wingate c 



ON 
OS 

< 



> 
o 

oo 

o 

On 



X 



Department of Physics, Washington University, St. Louis, MO 63130, USA 
b Department of Physics, Brookhaven National Lab, Upton, NY 11973, USA 
c Physics Department, University of Colorado, Boulder, CO 80309, USA 
d Physics Department, University of Utah, Salt Lake City, UT 84112, USA 
°Department of Physics, Indiana University, Bloomington, IN 47405, USA 
f SCRI, Florida State University, Tallahassee, FL 32306-4052, USA 
s Department of Physics, University of Arizona, Tucson, AZ 85721, USA 
''Department of Physics, University of California, Santa Barbara, CA 93106, USA 

We describe the current status of the MILC collaboration computation of fg , fg s , fr> , fn s and their ratios. 
Progress over the past year includes: better statistics and plateaus at j3 = 6.52 (quenched), j3 = 5.6 (Nf = 2) and 
f3 — 5.445 (Nf = 2), new runs with a wide range of dynamical quark masses at (3 — 5.5 (Nf = 2), an estimate 
of the systematic errors due to the chiral extrapolation, and an improved analysis which consistently takes into 
account both the correlations in the data at every stage and the systematic effects due to changing fitting ranges. 



Existing and planned experimental measure- 
ments of B-B and B s -B s mixing do not constrain 
the Cabibbo-Kobayashi-Maskawa matrix without 
knowledge of the heavy-light decay constants Jb 
and /b 3 and the corresponding £>-parameters. 
This fact has led to a major effort in the lattice 
community to calculate these quantities 0] . 

We have been computing heavy-light decay 
constants with Wilson fermions over the past 
three years. The goal of the current stage of 
the computation is twofold: 1) to extrapolate the 
quenched results to the continuum and estimate 
all systematic errors within the quenched approx- 
imation, and 2) to estimate errors due to quench- 
ing by comparing, at fixed lattice spacing, results 
from quenched and from Np = 2 dynamical lat- 
tices. 

The key ingredients in the calculation have 
been described in Table |l| gives the lattice 
parameters. Computations on the largest lattices 
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have been performed on the 512-node and 1024- 
node Intel Paragon computers at Oak Ridge Na- 
tional Laboratory; Paragons at Indiana Univer- 
sity and at the San Diego Supercomputer Center 
have been used for smaller lattices. 

We use the hopping parameter expansion of 
Henty and Kenway S for the heavy quarks. The 
static-light decay constants then come "for free" 
from the computation. However, as explained in 
0] , the procedure is optimized for the heavy- light 
case, and the static-light results are not usable 
when the physical volume is too large (runs B, E, 
and L-P). A dedicated computation of the static- 
light decay constants for these lattices is now in 
progress ||. We use a multi-source technique, 
with relative wavefunctions taken from the re- 
sults of the Kentucky group ||. Since none of 
the weak-coupling quenched lattices (and only a 
subset of the coarser ones) are affected, the new 
running is unlikely to have much impact on the 
extrapolation of the quenched results to the con- 
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Table 1 

Lattice parameters. Runs F, G, and L-P use 
variable-mass Wilson valence quarks and two fla- 
vors of fixed-mass staggered dynamical fcrmions; 
all other runs use quenched Wilson quarks. 



name 


(3 (am q ) 


size 


# configs. 


A 


5.7 


8 3 x 48 


200 


B 


5.7 


16 3 x 48 


100 


E 


5.85 


12 3 x 48 


100 


C 


6.0 


16 3 x 48 


100 


D 


6.3 


24 3 x 80 


100 


H 


6.52 


32 3 x 100 


60 


F 


5.7 (0.01) 


16 3 x 32 


49 


G 


5.6 (0.01) 


16 3 x 32 


200 


L 


5.445 (0.025) 


16 3 x 48 


100 


N 


5.5 (0.1) 


24 3 x 64 


100 





5.5 (0.05) 


24 3 x 64 


100 


M 


5.5 (0.025) 


20 3 x 64 


100 


P 


5.5 (0.0125) 


20 3 x 64 


100 



tinuum. However, the Np = 2 results at stronger 
coupling may be significantly altered. At present, 
we thus base our estimate of quenching effects on 
a comparison with the results from run G, where 
the current static-light technique should be reli- 
able and indeed gives decay constants that are 
consistent with the preliminary results from Q . 

A major improvement in the computation over 
the past year is the consistent use of covariant 
fits at every stage of the data analysis. Previ- 
ously, only the fits of /q 9 a/ Mq q vs. l/MQ q took 
correlations in the data into account. As is well 
known Q, the presence of small, poorly deter- 
mined eigenvalues in the covariance matrix can 
make covariant fits very unstable. We therefore 
had been forced to use noncovariant fits to the 
raw correlation functions and in the chiral extrap- 
olations. 

Our new technique has many similarities to 
those proposed in ||, but has some advantages, 
at least in the current project. It is based on 
a standard approach in factor analysis 0. We 
first compute the correlation matrix (the covari- 
ance matrix normalized with l's along the diag- 
onal) and find its eigenvalues and eigenvectors. 
We then reconstruct the correlation matrix from 



the eigenvectors, but omitting those correspond- 
ing to "small" eigenvalues. The resulting matrix 
is of course singular. It is made into an accept- 
able correlation matrix by restoring the l's along 
the diagonal. Finally, the corresponding covari- 
ance matrix is constructed and inverted, and is 
used in the standard way for making the fits. 

The above technique interpolates smoothly be- 
tween ordinary covariant fits, where no eigenval- 
ues are omitted, and noncovariant fits, where all 
eigenvalues are omitted. Furthermore, because 
the correlation (as opposed to covariance) ma- 
trix is used, the eigenvalues are normalized, with 
the average eigenvalue always equal to 1. This 
allows us to make a uniform determination of 
which eigenvalues to keep, which is very impor- 
tant since we are dealing with thousands of fits, 
and it is impossible to examine each fit by hand. 
Our standard procedure is to keep all eigenval- 
ues greater than 1 ; this accounts for typically 90- 
95% of the total covariance. Indeed, when one 
changes how the covariance matrix is computed, 
(for example, by changing the number of configu- 
rations eliminated in the jackknife), the eigenval- 
ues smaller than 1 generally change drastically. 
The approach eliminates unstable, "pathological" 
fits completely. We check that the final results are 
not significantly affected when one keeps several 
more (or several fewer) eigenvalues throughout. 

The use of covariant fits, with a reliable 
X 2 statistic, has led to three other important 
changes in the analysis over the past year. First, 
we have found that the ranges in Euclidean 
time over which the smeared-local and espe- 
cially the smeared-smeared correlation functions 
were previously fit were rather optimistic, giving 
X 2 /dof ~ 1.5 with 40 or 50 degrees of freedom. 
We now go out further in time to get reasonable 
confidence levels; this has increased the statistical 
errors somewhat and also has tended to raise the 
central values by about 1 (old) sigma. 

A second change involves the chiral extrapola- 
tions of to 2 and j~7r • Fig. |l| shows a typical chiral 
extrapolation of /„.. The standard linear fit in 
1/k, although apparently reasonable, has a poor 
confidence level once correlations in the data have 
been taken into account. This feature also ap- 
pears in fits of to 2 vs. and has been noted 
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previously in the literature ||. Although there 
are no chiral logs in f n in quenched chiral per- 
turbation theory at one loop ||, nonlinear terms 
can enter at higher order. However, the nonlin- 
ear terms in Fig. [I] could just as easily come from 
0(a) effects. For example, one can generally get 
acceptable confidence levels for linear fits to / w 
simply by using a different quark mass definition 
(e.g., mi(ft) = ln(l + 1/2k — 1/2k c )) on the hori- 
zontal axis. (Unfortunately, the curvature in 
increases with this definition.) With just three 
light quark masses, we are unable to investigate 
such higher order effects in detail. We therefore 
choose linear fits in finding the central values of 
decay constants, and take the difference with the 
quadratic fits as a systematic error estimate. This 
is a significant (~ 8%) error which we did not in- 
clude previously. We note in passing that linear 
fits for heavy-light masses and decay constants as 
a function of light quark mass are generally ac- 
ceptable. 
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Figure f. Fits to f v vs. 1/k for run D. The lin- 
ear fit (solid line) has confidence level 0.03; the 
quadratic "fit" (dashed line) has no degrees of 
freedom. The "bursts" show the extrapolation 
to k c , which in turn is determined by the corre- 
sponding fit (linear or quadratic) to vs. 1/k. 

A third change appears in the interpolation of 
fq q yl Mq q in inverse meson mass to vtlb- A typ- 
ical plot is shown in Fig. ||. When all the fits at 
earlier stages of the analysis are performed covari- 
antly, the fits to the "heavier heavies" plus static 



generally have higher confidence levels than fits to 
the "lighter heavies" plus static and are therefore 
used for the central values for quantities involv- 
ing the b quark. On the other hand, when simple 
noncovariant fits are used in the earlier analysis 
stages, the two fits in Fig. || typically have com- 
parable (and somewhat lower) confidence levels. 
This was the case last year, and we then chose the 
"lighter heavies" fit for the central values. The ef- 
fect on fs of the new choice varies with f3 in the 
range -3 to +10 MeV. 




Figure 2. f Qq {M Qq )i vs. l/M Qq for lattice D. 
The solid line is a quadratic fit (conf. level = 0.87) 
to the octagons ("heavier heavies" + static); the 
dashed line is a quadratic fit (conf. level = 0.37) 
to the crosses ("lighter heavies" + static). 

The extrapolation of quenched Jb to the con- 
tinuum is shown in Fig. 0. Errors on the points 
are larger than before E], mainly because they 
now include, in addition to statistical errors, the 
effects of varying the fitting ranges of the raw cor- 
relators. We take the linear fit to all the quenched 
data (solid line) to give the central value; the 
dashed line gives an error estimate. An extrap- 
olation of the Np = 2 data is not yet feasible, 
but may become so when we complete 1) the new 
static computation on the coarser lattices (L-P) 
and 2) new dynamical fermion runs at (5 = 5.6, 
24 3 x 64. 

Figure || is a similar plot for fs s ■ Note that, 
although the difference at a ~ 0.5 GeV between 
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Figure 3. fg vs. a. Octagons are quenched data; 
crosses, Nf = 2. The solid line is a linear fit to all 
quenched points (conf. lcvcl= 0.66); the dashed 
line is a a constant fit to the three quenched points 
with a < 0.5 GeV (conf. level= 0.76). The ex- 
trapolated values at a = are indicated by bursts. 
The scale is set by / w = 132 MeV throughout. 

quenched and Np = 2 results for fg is not so 
clear, the corresponding difference seems reason- 
ably convincing for fg s . We are hopeful that the 
improvements mentioned in the previous para- 
graph will sharpen these differences (if indeed 
they are present). 

The analysis of systematic errors is largely un- 
changed from Lattice 95 (|]. As described above, 
however, the errors due to chiral extrapolation 
are now estimated and included. 

The (still preliminary) results are: 

f B = 166(11)(28)(14) f D = 196(9)(14)(8) 
Sb. = 181(10)(36)(18) f D . = 211(7)(25)(11) 

^ = 1.10(2)(5)(8) If*- = 1.09(2)(5)(5) 
lB Id 

where the first error includes statistical errors and 
systematic effects of changing fitting ranges; the 
second, other errors within the quenched approx- 
imation; the third, an estimate of the quenching 
error. Decay constants are in MeV. The estimate 
of the quenching error is crude at present and 
may well increase significantly when a continuum 
extrapolation of the dynamical fermion results is 
possible. 
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Figure 4. fg 3 vs. a. Points and fits as in Fig. ||. 
The solid line has conf. level= 0.37; the dashed 
line, conf. lcvcl= 0.48. 
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